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ABSTRACT 


The  properties  of  a  classical  magnetic  dipole  are  investigated.  The 
vector  potential  A  is  derived  ab  initio  and  from  this  the  explicit  form  of  the 
magnetic  induction  b  is  deduced.  It  is  verified  that  this  dipole  magnetic 
field  is  both  solencidal  (V .  B  =  o)  and  irrotational  (V  x  B  =  o)  so  that 
Maxwell's  equations  are  satisfied  for  steady-state  conditions.  These 
properties  also  lead  to  the  existence  of  a  scalar  potential  <t>  which  is 
explicity  derived. 
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PROPERTIES  OF  A  MAGNETIC  DIPOLE 


INTRODUCTION 


Magnetic  phenomena  differ  from  electric  phenomena  since  there  are  no 
free  magnetic  charges  cr  poles.  The  basic  entity  in  magnetic  phenomena  is 
thus  the  magnetic  dipole  -  the  magnetic  ar.alogue  of  the  electric  dipole 
consisting  of  two  opposite  charges  placed  very  close  together. 


Bulk  magnetic  phenomena  such  as  magnetization  are  generally 
described  in  terms  of  accumulations  of  magnetic  dipoles.  However  standard 
works  on  electromagnetism  (for  example  ref3.  I  1-3 !)' usually  describe  only  the 
simple  features  of  the  dipole  without  verifying  ail  its  properties.  This 
note  explores  the  physical  and  mathematical  properties  of  the  classical 
magnetic  dipole.  These  properties  enable  the  dipole  to  be  considered  as  one 
of  the  fundamental  magnetostatic  units. 


2.  THE  VECTOR  POTENTIAL  A 


The  magnetic  dipole  can  be  described  by  an  infinitesimally  small 
current  loop  such  as  shown  in  Fig.  l.  To  calculate  the  magnetic  induction  B 
due  to  this  current  loop  it  is  customary  to  first  calculate  the  vector 
potential  A  defined  so  that 


B  =  V  x  A 


(  1  ) 


The  loop  of  radius  a  is  oriented  in  tha  x-y  plane  with  its  centre 
at  the  origin  of  the  rectangular  cartesian  coordinate  system  and  the  vector 
potential  is  to  be  determined  at  a  point  P.  For  convenience  a  spherical 
polar  coordinate  system  is  also  used  with  its  origin  coincident  with  the 
cartesian  coordinate  system  (the  point  0  in  Fig.l).  The  following  derivation 
of  A  is  similar  to  that  given  in  standard  texts  (1-21  but  is  included  here  for 
completeness.  Due  to  the  freedom  of  gauge  transformations  for  the  vector 


l 


potential,  7. A  can  be  defined  as  zero  (the  so-called  Coulomb  gauge 
condition) .  In  this  case  the  vector  potential  A  at  point  P  is  given  by  the 
volume  integral 


A  =  n  /  4  it  [  j(  r '  )  /  (  I  r-r  I  )  dV' 
o  J 


<2> 


where  J  is  the  current  density,  r  is  the  radial  vector  defining  the  point  P 
and  r'  is  the  vector  defining  the  points  on  the  current  loop  over  which  the 
integration  is  to  be  carried  out.  The  current  density  J  has  a  component  only 
in  the  6-di recti  on 


J  =  14 ( cose ' ) 4 ( r ' -a ) /a  (3) 


and  hence  so  has  the  vector  potential.  With  no  loss  of  generality  and  to 
simplify  the  calculations  the  point  P  is  chosen  so  that  «  =  0.  Then  only  the 
rectangular  component  J  =  J  is  retained.  Now  A  is  given  by 

y  0 

A  (r,o)  =  il  l/4#a  f  r'2  dr '  dQcos* '  4  (  cos  e  '  )  4  (  r '  -a  )  /  (  I  r-r '  I  )  (4) 

4  o  J 

where  d£2'  =  smo'do'd*'.  The  integration  is  simplified  by  the  two  delta 
functions:  the  current  density  has  non-zero  value  only  on  the  loop  where 
cosa'  =  0  and  r'  =  a.  When  the  integration  is  carried  out  over  the  two  delta 
functions  the  vector  potential  reduces  to 


A  (  r ,  0 ) 

* 


f2w 

l  0 


cos<fr’d»' 


(5) 


For  points  P  where  r>>a  the  denominator  may  be  expanded  in  a  Taylor 

series  as: 


l  r-r-  I 


-1 


1/r  (1  +  (a/r)sinecos«- 


+  0(a2/r2) 


(  6  ) 


If  this  expression  is  substituted  into  equation  (5)  and  terms  of 
second  order  and  higher  in  (a/r)  neglected,  then  the  vector  potential  is  given 
by 


A  (  r ,  0 ) 

$ 


* 2  „ 

f  cos*'(l+  ( a/r ) sinocos* ' 1 d* ' 
J  0 


4»r 


(7) 


m  la 

0  f  2w  2  ,  . 

- —  sinocos  fd*  , 

2  J  o 

4»r 


(  8  ) 


2 


s  i  n  o , 


r  2  *  2 

since  ccs 
J  o 


o'd*'  =  *  and  f2*  cos«'d$'  =  0. 
J  o 


Defining  the  dipole  moment  m  to  be  a  vector  of  magnitude 

( ia  )l  normal  to  the  plane  of  the  current  loop,  A  <r,o)  can  be  written  as 

$ 

,  ,  'o 

A  ( r , o )  =  —  msino/r 
0  4  * 

Hence  the  vector  potential  A  can  be  expressed  as  a  vector  cross 

oroduct 


A  -  n  / ir  m  x  r/r 3 , 

0 


since  its  only  non-zero  component  (in  the  ^-direction)  is  orthogonal  to  both 
the  directions  of  m  and  r  and  its  magnitude  is  given  by  ?  /4*  imMr/r3! 
times  the  sine  of  the  angle  e  between  these  two  vectors. 

Alternatively  A  can  be  expressed  as: 

A=-(i/4*mx7(l/r)  (12) 

c 

since  7(l/r)  =  =  -r/r3  for  r  >  0. 

If  this  expression  for  A  is  correct  then  it  should  satisfy  the 
Coulomb  gauge  condition  assumed  in  its  derivation.  This  follows  readily  since 


7. A  =  it  / 4 *  7  •  ( m  x  r/r' 


/ 4 *  ((7  x  m)  r/r3  +  m.(v  x  r/r3))  =  o 


because  7  x  m  vanishes  as  m  is  constant  and  the  curl  of  the  gradient  of  any 
function  vanishes  identically  (since  r/r3  =  -  7 ( 1 / r ) > . 


3 


3.  COMPONENTS  OF  MAGNETIC  INDUCTION  B 

T'r.e  magnetic  induction  B  due  to  a  magnetic  dipole  located  at  the 
origin  can  be  derived  from,  the  relationship 

B  =  V  x  A  (15) 


If  the  result  obtained  in  section  2  for  the  vector  potential  A  is 
substituted,  the  magnetic  induction  at  a  point  r  away  from  the  current  source 
is  given  by: 

B  =  ,t/4ir7x(mxr/r3) 


which  can  be  expanded  to 

B  =  n  /  4  ir  [(r/r3.7)m  -  (?.m)r/r3  +m(7.r/r3)  -  (m.v)r/r3]  !  1 ~  > 


The  first  two  terms  are  identically  zero  since  m  is  constant.  The  third  term 
is  also  zero  since  for  any  integer  n  and  for  r  >  0  (4) 

7 . ( rnr)  =  rn ( n  +  3  >  ! :  3 ! 


Thus  the  only  remaining  term  is 

B  =  -  p  /4»  (m. 7) r/r3  (19) 

o 

This  expression  can  be  simplified  by  considering  the  vector 

I  ipr.!!!  ty ! 

(0-7)f  V  =  (UV^)V  +  *(U'V)V  (20) 


where  U  and  V  are  arbitrary  differentiable  vector  functions  and  #  is  an 
arbitrary  differentiable  scalar  function.  Equation  (19)  then  becomes 


B  =  — -  [ ( m- 7 ( l /r 3 ) ) r  +  - —  (mv)r) 
4  *  3 

K 

~  M 

*  7 —  [-3(m  r)r/r5  +  m/r3l 

4 » 


(21) 


!  22) 


4 


because  v(i/r3)  =  -  3r/r5  and  (m-V>r  =  m.  Thus  (22)  can  be  written 


B  =  —  [3(m-r)r/r5  -  m/ r3  . 
4  IT 


(23) 


In  spherical  polar  coordinates  with  the  origin  situated  at  the  location  of  the 
dipole  the  components  of  B  can  be  found  from  (23)  by  Br  =  r •  B  and  B0  =  o  •  B 
where  r  and  o  are  unit  vectors  in  the  radial  and  zenithal  (o)  directions 
respectively  (ref.  to  Fig.  l).  There  is  no  component  in  the 
azimuthal  (*)  direction  since  «■  B  =  0.  The  two  non-zero  components  may  be 
readily  obtained  as: 


3  =  n  /  4  ir  !2mcoso/r)  (24) 

r  o 


3  =  f.  / 4 *  (msino/rl  (25) 

8  o 

from  the  above  discussion. 


4.  DIVERGENCE  OF  THE  MAGNETIC  INDUCTION 


The  expression  (23)  for  B  was  derived  from  the  curl  of  the  vector 
potential  A.  Hence  b  is  automatically  divergenceless  since  the  divergence  of 
the  curl  of  any  vector  is  zero  14!.  However,  it  is  of  interest  to  explicity 
verify  that  the  deduced  expression  for  B  is  divergenceless.  The  divergence  of 
the  magnetic  induction  due  to  a  magnetic  dipole  is  given  by 


7  B  =  —  7 ■ ( 3 ( m ■ r)r/r5 

4  r 


m/r  3  ) 


(  26  ) 


The  divergence  of  both  terms  inside  the  bracket  can  be  calculated 
from  the  vector  identity  (4) 


7 ■ ( fv)  =  V*.V  +  *VV  ( 27) 

where  f  and  V  are  defined  as  in  section  3.  Now  the  first  term  may  be  expanded 
to 

7.(m.r)r/r5  =  m.r/r5  +  ( m. r ) ( 7 . r/r5 )  (28) 

=  -m. r/r 5  (29) 


5 


c;  5 

since  7(m.r)  =  m  and  7. (r/r")  =  -2/r  from  equation  ( 1 8 ) .  The  divergence  of 
the  second  term  in  (26)  reduces  to 


7.(m/r3)  =  V(l/r3).m  +  l/r3(7.m) 

/  5 

=  -3m. r/r 


since  v(i/r  )  =  -3r/r  and  m  is  a  constant.  Combining  these  results  the 
divergence  of  B  becomes 

«  5  5 

7 . B  =  —  !-3m.r/r  +  3m.r/r  1  <32/ 

4  T 


=  0 , 


so  that  B  is  a  sclenoidal  vector  field  and  thus  obeys  the  fundamental  law  of 
magnetostatics . 


5.  THE  CURL  OF  THE  MAGNETIC  INDUCTION 


From  Maxwell's  laws  v  x  B  =  n  J  assuming  steady-state  conditions. 
Hence  v  x  B  =  0  in  the  region  excluding  the  dipole  source.  Adapting  the  same 
philosophy  as  in  section  4  this  result  is  explicity  verified  for  the  specific 
form  of  the  dipole's  magnetic  induction  B.  The  curl  of  the  magnetic 
induction  can  be  written  as 


7  X 


—  7  x  (3<m.r>  r/r5 

4  K 


m/r 3  ) 


(  33  ) 


The  curl  of  both  terms  can  be  calculated  from  the  vector  identity 

7  X  ( *V)  =  7#  X  V  +  f(7  X  V),  (34) 


where  f  and  V  are  as  defined  in  section  3.  The  first  term  may  be  expanded  as 


37  X 


( m.  r)  r/  r 


5 


5  5 

31  7(m.r)  x  r/r  +  (m.r)  7  x(r/r  >) 

=  3 1m  x  r/r5 1 


(  35  ) 

(  36 ) 
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The  second  term  m 


since  7(m.r)  =  m  and  7  x  (r/'r5)  =  0  identically  Ml. 
equation  (33)  becorr.es 


-7 


x  ( m/rJ ) 


3  3 

[7(i / r  )  xm+  l/r  (7xm 


(37) 


=  3 


r  x  m/r 


5 


l  38) 


3  5 

since  7(l/r  )  =  -3r/r  for  r  >  0  and  v  x  m  =  0.  combining  these  results  the 
curl  of  B  becomes 


7  x  B 


[  3 (m  x  r) /r 5 


3  (  r  x  m) /r 5 ] 


(  39 ) 


Hence  the  magnetic  induction  vector  B  is  irrotaticnal  as  required. 


6.  SCALAR  POTENTIAL  OF  MAGNETIC  DIPOLE 


It  has  been  shown  that  B  is  irrotational  (i.e.  7  x  B  =  0)  for  a 
magnetic  dipole.  Hence  B  can  be  written  as  the  gradient  of  a  scalar 
field  <t> . 


B  *  -  7* 


(  40  ) 


Since  g  is  also  solenoidal  (V  B  =  0)  the  scalar  potential  satisfies 
the  Laplacian,  7%  =  0. 

The  scalar  potential  can  be  deduced  from  the  alternative  expression 
for  the  vector  potential  A 


A  =  -*  /4»  m  x  v( l/r)  (41) 

0 

Thus  B  =  7  x  A  can  be  expanded  as 


7 


B  =  —  tm.7)7(l/r)  -  m. ( 7 . 7( l/r ) ) 1  '42! 

4  ir 

=  —  [ (m.7 ) v( l/r )  ]  ,  (43) 

4  * 

2 

since  7  ( l / r )  =  0  for  r  >  0. 

New  using  the  vector  identity  [41 

7 ( m. 7 { l/r ) )  =  ( 7 ( l/r ) . 7 ) m  +  (m.7)v(l/r)  +7<l/r)  x  (7  x  m)  + 

mx(7x7(i/r>)  >'44! 


it  car.  be  seen  that 


7 ( m.7 ( l/r) )  =  ( m. 7 ) 7 ( l/r) 


!  45  ) 


since  the  other  terr.c  on  the  right  side  vanish  because  m  is  constant  and 
7  x  7(l/r)  vanishes  identically  !4l.  Thus  B  can  be  expressed  as 


B  =  —  7  [ (m.7< l/r) ! : 


so  that  the  scalar  potential  «  is  given  by 


«  =  -  «  /4t  m.7( l/r) 
o 


(47) 


It  can  be  readily  shown  that  this  potential  satisfies  the  Laplacian 


since 


7  «  =  -7  .  B  =  0 


7.  SUMMARY  AND  CONCLUSIONS 


In  this  paper  the  magnetostatic  properties  of  a  magnetic  dipole  have 
been  verified.  In  summary  it  has  been  shown  that  the  magnetic  dipole  has  a 
vector  potential  which  is  given  by 


8 


■  3-, 

<  m  x  r/  r  ) 


A  =  —  ;  m  x  7 ! ;  /  r  !  ) 


(50) 


ar.d  satisfies  the  Coulomb  gauge  cor.ditior 


7 .  A  =  o 


!  51  ) 


The  magnetic  induction  B  was  derived  from  the  curl  of  the  vector 
:er. tiai  to  be 


B  =  —  [  3  ( m .  r )  r/  r  5  -  m/  r  3 

4  * 


(  52  ! 


It  has  also  been  explicitly  verified  that  Maxwell's  equations  hold 
for  the  dipole  field  sc  that 


7  .  B  =  0 


(53/ 


an: 


7  X  B  =  0, 

the  region  excluding  the  current  source. 

Further,  it  was  shown  that  B  can  be  expressed  as  the  gradient 
:alar  potential  defined  by 


*  =  —  m.7  ( l/r ) 

4  » 


(  5  - 


;f  a 
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FIGURE  1 


The  magnetic  dipole  represented  as 
both  its  rectangular  cartesian  and 
systems . 


a  circular  current  loop  with 
spherical  polar  coordinate 
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